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Abstract 

We analyze the spectral properties of the volume operator of Ashtekar and Lewandowski in 
Loop Quantum Gravity, which is the quantum analogue of the classical volume expression for 
regions in three dimensional Riemannian space. Our analysis considers for the first time generic 
graph vertices of valence greater than four. Here we find that the geometry of the underlying 
vertex characterizes the spectral properties of the volume operator, in particular the presence 
of a 'volume gap' (a smallest non-zero eigenvalue in the spectrum) is found to depend on the 
vertex embedding. We compute the set of all non-spatially diffeomorphic non-coplanar vertex 
embeddings for vertices of valence 5-7, and argue that these sets can be used to label spatial 
diffeomorphism invariant states. We observe how gauge invariance connects vertex geometry 
and representation properties of the underlying gauge group in a natural way. Analytical results 
on the spectrum on 4-valent vertices are included, for which the presence of a volume gap is 
proved. This paper presents our main results; details are provided by a companion paper [52J. 

1 Introduction 

Loop Quantum Gravity (LQG) [H EJ [3j has become a promising candidate for a quantum theory of 
gravity over the last 15 years. It is an attempt to canonically quantize General Relativity (GR) while 
preserving its key principle: background independence^. The resulting quantum theory is formulated 
as an SU (2) gauge theory. 

* brunnemann® mat h . uni- hambur g . de 
I'd. rideout@imperial.ac.uk 

1 With the term 'background independence' we imply independence on a choice of fixed background geometry. In 
order to allow for the possibility of topology change in quantum gravity, it would be desirable to have the theory in 
its final formulation be independent of the topology of the underlying manifold, however to our knowledge this has 
not been achieved so far within LQG. 
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Upon casting GR into a Hamiltonian formulation by introducing a foliation of four dimensional 
spacetime (A^,(? M ^j§ into spatial three dimensional hypersurfaces E, with the orthogonal timelike 
direction parametrized by t e R, one obtains first class constraints which have to be imposed on the 
reformulated theory such that it obeys the dynamics of Einstein's equations and is independent of 
the particular choice of foliation. These constraints are the three spatial diffeomorphism or vector 
constraints which generate diffeomorphisms inside S, and the so called Hamilton constraint which 
generates deformations of the hypersurfaces E in the t- (foliation) direction. In addition one ob- 
tains three Gauss constraints due to the introduction of additional SU (2) gauge degrees of freedom. 
(There is also a more recent reformulation of the constraints, as proposed in the master constraint 
programme. See [9] and references therein.) The resulting theory is then quantized on the kinemati- 
cal level in terms of holonomies h and electric fluxes E. Kinematical states are defined over collection 
of edges of embedded graphs. The physical states have to be constructed by imposing the operator 
version of the constraints on the thus defined kinematical theory. 

There exist well defined operators in the kinematical quantum theory which correspond to (classi- 
cal) differential geometric objects, such as the length of curves [TO] , the area of surfaces [41 J, and the 
volume of regions [341 143] in the spatial foliation hypersurfaces E. Remarkably all these quantities 
have discrete spectra, which can be traced back to the compactness of the gauge group SU{2). 

The thus constructed kinematical theory has already proved to potentially answer a number of 
outstanding physical questions. Matter can be coupled to gravity without any problem, and the ul- 
traviolet divergences which plague conventional quantum field theory (QFT) on a fixed background 
disappear if geometry and matter are simultaneously quantized [8]. A coherent state framework 
[TT] has been formulated which allows one to address questions on the correct semiclassical limit of 
the theory as well as to examine the possible effect of a dynamical background on the propagation 
of matter [TJJ H3j- Moreover, the entropy of a black hole has been calculated from first principles 
PSl Q51 E] in the LQG framework. 

A central role in all these investigations is played by the area and the volume operator respectively. 
The former is understood very well analytically [4T] and numerically [47] and enters the construction 
of coherent states [121 [131 E] and the computation of black hole entropy [EEl EETJ- However, less 
is known about the latter, due to its complicated structure. Nevertheless it is a crucial object not 
only in order to analyze matter coupling to LQG, but also for evaluating the action of the Hamilton 
constraint (alternatively master constraint) operator in order to construct the physical sector of LQG. 
By construction the spectral properties of the constraint operators are driven by the spectrum of 
the volume operator. In addition, the analysis of the classical big bang singularity [2S1 EE] heavily 
depends on the volume operator. As is obvious, controlling this operator is highly desirable, if one 
wants to truly test the present framework of LQG. In particular it is of utmost importance to check 
whether the promising results obtained in symmetry reduced frameworks (see e.g. [191 1201 1221 123] 
and references therein) remain valid in full unreduced LQG. From this perspective the present work 
is also intended to serve as a starting point for numerical studies which explore physical predictions 
from the full theory. 

Historically two quantizations of the classical volume expression have been worked out: In [37] 
Loll has constructed a quantized volume expression on a fixed lattice. This expression was later 
slightly changed [38] and in [39] recognized to be compatible with the construction by Ashtekar 

2 Here denotes the four dimensional metric on M. which is a solution to Einstein's equations. Note that in this 
context, because of the diffeomorphism invariance of General Relativity, the metric is a representative of an equivalence 
class of metrics which can be transformed into each other by spacetime diffeomorphisms. 
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and Lewandowski y( AL ) in [42J and later in [33]. The second construction V( RS ^ [34] due to Rovelli 
and Smolin uses a different regularization scheme. These two operators differ in their invariance 
properties: while the operator V^ RS ^ is invariant under homeomorphisms of the spatial foliation 
hypersurfaces E, the version V^ AL > is invariant only under diffeomorphisms in E, its spectrum is 
sensitive to non-differentiable transformations. 

The action of V^ RS ^ on the kinematical basis states has been analyzed in [35], in [36] its matrix 
elements are computed for the special case of four-valent vertices and in [441 143] bounds on the 
eigenvalues for general four-valent vertices and monochromatic n-vertices are given. A closed formula 
for the general matrix elements of ~\f( AL \ as well as eigenvalues for special spin configurations, have 
been obtained in [40]. However, the derived expression suffers from its complicated structure which 
in practice prevents a more detailed analysis. The situation was significantly changed with the results 
of [S_TJ, in which a drastically simplified expression for the matrix elements of V^ AL ' was obtained. 
This opens the way to analyzing V^ AL ' in general as presented below. 

In [49l [50] the two operator versions V^ RS \ have been compared with respect to their 

consistency upon using an alternative electric flux quantization!!. The result is that using V^ AL ) is 
consistent, while the usage of leads to problems. For this reason we employ the regularization 

used in the construction of V^ AL -*, and imply V^ 1 ^ when referring to the volume operator in the 
remainder. 

Of course, one may ask several questions at this point: Why should we expect a quantum theory 
of gravity to be sensitive to the differentiable structure of the underlying (foliation) hypersurfaces? 
Shouldn't we regard this as a structure only to be re-obtained in a classical limit of the ultimate 
theory? Also, why should we consider vertex valences greater than four? One can introduce a 
triangulation of a three dimensional foliation hypersurface E by decomposing it into tetrahedrons. 
This cell decomposition is then dual to a graph having four valent vertices only, one vertex inside 
each tetrahedron. As the volume operator acts only non-trivially at the vertices (see below), we can 
well interpret this as a cell decomposition of E, where each cell comes with a volume according to the 
action of the volume operator. Moreover the Hamilton constraint operator as constructed in [H] does 
not increase the valence of the vertices of a graph, so the four-valent sector seems to be preserved 
under 'evolution' in foliation direction t. 

We would like to take a rather conservative point of view here: By construction LQG is a 
minimalist theory in the sense that, starting as close as possible to General Relativity, one introduces 
as few additional structures as necessary in order to construct a quantum theory of gravity. One then 
wishes to analyze the fate of classical quantities in the quantum theory. The quantization in terms of 
holonomies and fluxes seems to be a natural choice in order to quantize in a background independent 
way. We can construct a representation of the classical holonomy-flux-algebra on the kinematical 
Hilbert space, where holonomies act by multiplication, fluxes by differentiation, on states. Here it 
is important to note that multiplication is well defined for holonomies in an arbitrary graph, which 
contains arbitrary valent vertices. Our intent is to study the predictions of the theory as it is stated 
formally, without making any assumptions regarding the nature of physical states. 

Of course, the construction in terms of graphs suggests an interpretation of the kinematical the- 
ory in terms of a triangulation of E, hence it is tempting to directly associate graphs as parts of a 
triangulation. However, one also might take the point of view that graphs are only abstract auxiliary 

3 The idea is to regard the volume operator as fundamental and to derive an operator corresponding to the flux 
operator from the volume. One may then check whether the thus constructed alternative flux quantization is in 
agreement with the usual one. See section [5] 
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constructions one must introduce in order to describe a meaningful quantum theory, but without a 
direct interpretation. It would then be desirable to derive a possible interpretation from the abstract 
Hilbert space formulation rather than introduce it by hand. In particular one should be able to 
show from first principles of the abstract formulation that it is sufficient to restrict to certain (e.g. 
4-valent) types of graphs. 

With this motivation at hand we have numericalljQ analyzed the spectrum of the volume operator 
at non-coplanar vertices (no three edge tangents lie in a single plane) of valences 4-7, and also the 
widely used [13J, [18] cubical 6-vertex. We find that the particular embedding of a graph into E on 
the one hand and the recoupling of representations of the gauge group SU{2) on the other hand 
are combined in the spectral properties of the volume operator. The distribution of eigenvalues at a 
given vertex is driven by the embedding of this vertex, expressed in terms of (spatial) diffeomorphism 
invariant sign factors which label the relative orientations of edges at the vertex. In particular the 
existence of a volume gap depends on the vertex embedding. We find that eigenvalues accumulate 
close to zero in general. For the cubic six-valent vertex, however, we observe the presence of a volume 
gap. Also, for arbitrary gauge invariant four-valent vertices, we prove the existence of a volume gap 
and give an analytic lower bound for the spectrum. 

The diffeomorphism invariant sign factors, moreover, allow for a classification of possible non- 
diffeomorphic embeddings of a vertex having N edges, and can be used as a discrete label for diffeo- 
morphism invariant states^. We will show how the possible embeddings of an iV-valent vertex can 
be computed numerically and analytically. As it turns out, there is an important interplay between 
gauge invariance and sign factors which leads to cancellations in the volume contributions for certain 
large valence vertices. This suggests that it may be possible to formulate conditions for which higher 
valent vertices can indeed be neglected. A more detailed analysis will be presented in a forthcoming 
paper [53] . 

This paper is organized as follows. In section 2 we will recall basic steps in the construction of 
LQG as necessary in order to explain the starting point for our analysis. At the same time this will 
introduce our notation. Note that this introductory chapter is not reproduced in the companion 
paper [52]. Experts may safely skip this section. The construction of the volume operator V( AL ^ is 
sketched in section 3, in particular the occurrence of the sign- factors is illustrated there. We sketch 
the implementation of the volume operator on a computer and discuss the problem of computing 
gauge invariant basis states as well as the set of possible non-diffeomorphic embeddings of vertices 
into the three dimensional foliation hypersurfaces E. Section 4 provides the results of a detailed 
numerical spectral analysis of the volume operator for vertex valences 5, 6, 7, where the sensitivity 
of V to the diffeomorphism invariant geometric properties of the underlying vertex is demonstrated. 
Section 5 finally states analytic results on the volume operator action at four-valent vertices, in 
particular we give an analytic lower bound for the smallest non-zero eigenvalues. We close with a 

The structure of the matrix representation of the volume operator becomes very complicated already for vertices 
of valences of 5 and above. This in general prevents an analytical treatment. See however section [5] for analytical 
results on the 4-vertex. 

5 This vertex can be compared to the lattice points in a cubic lattice. Certainly opposite edges have parallel tangents 
in this case, hence there are triples of coplanar edge tangents. 

6 As shown in [27], there also exists a set of continuous diffeomorphism invariant parameters, called moduli, which 
additionally label a vertex. However the volume operator is not sensitive to this information, and classifying states 
with respect to their volume spectrum results in the described discrete set of vertex classifications labelled by sign 
factors. 
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summary and outlook section 6 where we point to future application of our results, as well as the 
computational tools developed during this project. 

We close the introduction with a side remark. The present publication summarizes our most 
important results. The interested reader is referred to [52] to study the details on each of the results 
presented here. We recommend reading this summary paper first in order to get an overview of what 
has been achieved, before studying the lengthy detailed presentation in [52J. 

2 Loop Quantum Gravity 

2.1 Classical Theory 
2.1.1 Hamiltonian Formulation 

In order to cast General Relativity into the Hamiltonian formalism, one has to perform a foliation 
M = R x S of the four dimensional spacetime manifold (M.^^) into three dimensional spacelike 
hypersurfaces S with transverse time direction labelled by a foliation parameter t e R. Note that 
this foliation has to be seen as an auxiliary construction and it is kept completely arbitrary, since 
General Relativity is manifestly invariant under four dimensional diffeomorphisms of spacetime. 

As analyzed by Arnowitt, Deser and Misner [3], the four dimensional metric g^ u can then be 
decomposed into the three metric q a b(x) on the spatial slices 10 and its extrinsic curvature -Kafc(z/jE 
which serve as canonical variables, (x and y are points in £.) However it seems to be difficult to 
formulate a quantum theory out of this approach. These problems can be overcome if one introduces 
new variables due to Ashtekar [5] : one may implement additional gauge degrees of freedom and take 
as canonical variables electric fields Ef(x) and connections A 3 b (y) as used in the canonical formulation 
of Yang Mills theories. Here a,b — 1 ... 3 denote spatial (tensor) indices, i,j = 1 . . . 3 denote su(2)- 
indices. 

The occurrence of SU(2) as a gauge group comes from the fact that it is necessary for the coupling 
of spinorial matter and it is the universal covering group of 50(3) which arises naturally when one 
rewrites the three metric q a b in terms of cotriads as 

Qab = e{e k b 5 jk . (2.1) 

Obviously this definition is only unique up to 50(3) rotations. 
The pair (A, E) is related to (q, K) as 

A{ = rj - pK ab e) E« = V / dit^)e« (2.2) 

with T being the spin connection (a certain function of the triads e{ and its first spatial derivatives) 
and e| is the inverse of e J a , that is e°-e\ = S b , e"e^ = 5j. Moreover we have the so called Immirzi 
parameter (3 [Bj, which in principle can be chosen as an arbitrary complex constant not specified a 

7 As induced by g^, that is the projection of g^ v on the spatial slice. 

8 That is the rate of change of q a b{x) with respect to t. This encapsulates the information of g^ v in the foliation 
direction. 
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priori, but is chosen to be real in most caseg^|. The pair (A, E) obeys the Poisson bracket: 

{Et(x),Ai(y)} = K 5 a b 5{5(x,y) (2.3) 

with k = 8ttGn, Gn being Newton's constant. All other Poisson brackets vanish. Note that (12.31) is 
invariant under a rescaling by the Immirzi parameter as Ef(x) — ► jjE°-(x), A b {y) — ► f3A 3 b (y). For 
clarity of the notation we will suppress this j3 dependence in the following and only insert it where 
necessary. 



2.1.2 Constraints 

The constructed theory is subject to constraints which arise due to the background independence of 
General Relativity and the resulting non-bijective nature of the Legendre transformation performed 
upon introducing conjugated momenta. It can be treated as suggested by Dirac [57]. One obtains 
the three vector (spatial diffeomorphism) constraints 

Ca = F 3 ab E) (2.4) 

where F 3 ab = d a A 3 h — dbA 3 a + ejkiA^A l b is the curvature of the connection A. Moreover one gets a 
scalar(Hamilton) constraint 



C — — = £jki 



FL - (1 + P 2 )e 3 mnKK% E a k E\ (2.5) 



b 



with K{ — A J a — F 3 a (E). Additionally, due to the introduction of additional SU{2) degrees of freedom 
in ( 12. ip . one obtains the three Gauss constraints 

Cj = d a E] + e jkl A k a E? (2.6) 



2.1.3 Regularization of the Poisson Algebra 

We have shown how General Relativity can be rewritten as a theory with constraints in the Hamilto- 
nian formalism. Moreover we have chosen a polarization of the classical phase space, that is we have 
chosen a split into a configuration space A, consisting of smooth S'L r (2)-connections A J b (x) defined on 
every point x of the spatial foliation hypersurfaces n 10 i and a momentum space T*A, the cotangent 
bundle containing the momenta Ef(x). 

As mentioned above we must enlarge the phase space by introducing SU{2) gauge degrees of 
freedom, that is we can perform gauge local SU(2) transformations g 6 Q, where Q denotes the 
space of smooth gauge transformations and every g is defined locally as g(x) for all x G S. 

Now the Poisson algebra (12. 3p has to be regularized in order to construct a meaningful quantum 
theory. Therefore one uses the following smearec0 quantities : 

9 It is possible to set (3 to a complex number, which tremendously simplifies the form of the constraints below. 
However this simplification comes at the price of ensuring the reality of geometric quantities such as the spatial metric 
q a b- It turns out that these reality conditions are very difficult to implement in the quantum theory. This is the reason 
why [3 is usually chosen to be real. 

10 In the following we will refer to an A 6 A as a particular smooth configuration, given by its values Al(x) at every 
x € E. 

11 For details on the regularization procedure we refer to [Tj. 
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The duao *Ej(x) (a 2-form) of Ej(x) integrated over a two dimensional orientable piecewise 

analytic surface S, embedded in E by 

(S)X ■ [ ~ 1 \] 2 3 ( u \ u 2 )^S = X{u\ u 2 ) : 

E t (S) := J*E l = J d 2 u e a J s X b ul {u\u 2 ) ^X% 2 {u\u 2 )E-( S) X{u\u 2 )) (2.7) 



The integrated version of A J b (y)(a, 1-form) over a one dimensional oriented piecewise analyti 
edges e embedded in £ by 

e : [0, 1] 3 1 1— > e(t), that is the holonomy h e (A) of A 3 h {y) given by 
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i i i 

oo 



h e (A) :=Vexp J A =1 2 + jr J dh J dt 2 ... J dt n A(ti)A(t 2 ) . . . A(t n ) (2.8) 

n ~ 1 ti tn-l 

Here Pexp- denotes the path ordered exponential. Both Ei(S) and h e (A) are by construction 
invariant under reparametrizations of the surface S C S and the edge ecS. When we use the term 
'edge' in the following we will actually imply an equivalence class of edges, which differ only by a 
reparametrization. Such an equivalence class is sometimes also referred to as a path. 

Note that an arbitrary element g of the set Q of smooth gauge transformations gives a map 
g : T, 3 x ^ g(x) G SU(2), that is to every point x G S an element g(x) of the gauge group SU(2) 
is assigned. As the connection A transforms under g as A — > A 9 := ad g A — dgg" 1 the holonomy h e 
transforms under g as : 

hi{A) := h e (A 9 ) = 9{b{e)) K{A) ^(/(e))" 1 (2.9) 

with 6(e) = e(t = 0), /(e) = e(t = 1). Now assume that the edge e intersects S at the intersection 
point e(ti nt ) given by the curve parameter t int . For an arbitrary edge e one can always obtain a 
splitting e = ei o e 2 such that the edge e is decomposed into two edges ei, e 2 . Accordingly we then 
have h e = h ei ■ h e2 where • denotes group multiplication. Moreover we can change the orientation of 
ei,e 2 such that their common starting point is e(t int ), then h e -i = (Zie)" 1 . Therefore without loss of 
generality one can restrict oneself to consider only edges outgoing from e(t int ). Then (12. 3p can be 
rewritten as: 

{Ei{S),h e (A)} = | e(e,S) r 3 h e (2.10) 

Here the sign factor e(e, S) = if e, S do not intersect or if e C S. If e, S intersect in a point given 
by the curve parameter t int then e(e, S) = 1 if the normal of S and the tangent vector e(ti nt ) point 
to the same side of S and e(e, S) = —1 if they point to different sides. In what follows we always 
will assume such an adaption of edges to the surface S when acting with Ei(S). 



2.1.4 Cylindrical Functions 

Let us now consider a graph 7 embedded in the spatial hypersurface E, that is a collection £'(7) : = 
{ei, . . . , e^} of I £(7) I := iV edges (ei, . . . , ejv), having common points at their endpoints only. (N 

12 With respect to spatial degrees of freedom. 
13 Edges should intersect only finitely many times. 
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can be countably infinite.) The set of these intersection points of edges is called the vertex set V(^y) 
of the graph 7. Moreover let us denote by T the set of all such graphs 7 C S. 

The set of holonomies {/i e } eg £( 7 ) then defines a map p 7 from the space A of smooth connections 
A to the N copies of the gauge group SU (2) via: 

p 7 : A -> Sf/(2)l^l 

A i-> {h e (A)} eeEM (2.11) 

A function / : A — > C is then called cylindrical with respect to a graph 7 if it can be written as 
/ = fj o p 7 such that 

/ 7 (MA),...,M^)) : 5C/(2)I E «I^C (2.12) 
where / 7 depends on the holonomies of the edges e G E{p(), each defining a group element on one 
of the N copies^ of the group SU(2). By applying Leibnitz' rule in the evaluation of the Poisson 
bracket (12.101) one obtains: 

{E l (S)J 1 ({h e (A)}eeE b) )} = ^ £ E e ( e ' 5 ) [nhe] AB g^T-M{he(A)} eeE{l) ) (2.13) 

Note that jr^/ie] AB d(h 9 ) A B can ^ e identifiecl with the action of the right invariant vector field (Xjf) (g) := 
■^f(e tTj g)\ t =o on G, that is, the vector field constructed by the transport of the tangent space T t G 
at the unit element 1q of G to all elements g G G via right-multiplication with g. 



2.2 Quantum Theory 

2.2.1 Flux and Holonomy Operators 

Now formally the Poisson bracket {Ej(S), /i e } can be turned into ih[Ej(S), h e ] , the commutator of 
the according operators multiplied by ih. This gives us Ej(S) acting as a differentiation and h e being 
the multiplication operator on the space of cylindrical functions CYL of the form / = / 7 o p 7 . We 
will often suppress the decorations on h e . 

Note that the formulation in terms of holonomies (12. 8p and fluxes allows for connections A G A 
contained in the space of distributional connections A. In other words: Due to the compactness of 
SU{2) the cylindrical functions are bounded and hence to each cylindrical function we can associate 
its sup-norm. Then the closure CYL with respect to the sup-norm can be associated to A, which can 
be seen as the set of arbitrarily discontinuous maps from the set of edges (paths) to the gauge group 
SU(2). Moreover the gauge behavior (I2.9p allows for gauge transformations g G Q contained in the 
set of distributional gauge transformations Q. To summarize: A l a (x) and g(x) may vary arbitrarily 
within an infinitesimal neighborhood U(x) of a point x G Moreover we will denote by A\^ C A 

the restriction of the connection configurations A G A to the one dimensional subset 7 C E. Then 
Q23SD gives: 

EJ(S) [/ 7 ] = ih{E t (S)J\({h e (A)UE { ,))} =J4E <^ S ) X l A ( 2 - 14 ) 

14 One copy for each e e E(j). 

15 In contrast to the classical theory which considers only smooth and even differentiable quantities which can be 
described by their Taylor series, this means that we have assignments A 9 A : x — > A l a (x) and Q 3 g : x — > g(x) for 
every point x in S, with no further restrictions such as differentiability or smoothness. 
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where we have denoted the right invariant vector field on the copy of SU(2) belonging to a particular 
edge e G E{^y) by X % e1 again % = 1,2,3 being an su(2) index 16 !. Here naturally the Planck area 
shows up by £ 2 P = Hk. Furthermore we pinned an index 7 to E to denote that this operator 
acts on a function / 7 restricted to the graph 7. One is then led to consider the Hilbert space 
H y := L 2 (A\ ,dfM~f). 7i 7 is equipped with a natural inner product constructed from the measure 
dfi-f := riee-B(7) dnH(h e ) (the product of the Haar measures d\iu on the individual copies of SU(2) 
on the edges e G E ('-/)). It contains (the completion with respect to the sup-norm of) the set of 
cylindrical functions / 7 (/i ei (v4), . . . , h eN (A)), A G A\ , square integrable with respect to dfi^. 

2.2.2 Choosing a Basis in 7i 7 : Spin Network Functions 

Recall that irreducible representations of SU{2) can be realized by unitary square matrices [^-(/i)] 
acting on a (2j + l)-dimensional linear vector space. The label j is a non-negative half integer 
(j = 0, |, 1, . . .) and is called the weight of the representation ttj. Moreover matrix elements [^-(/i)] mn 
of [vrj(/i)] are usually labelled by m,n — j, ... , —j from the upper left to the lower right corner. 

Now due to the Peter & Weyl theorem the matrix element functions [tTjO)] °f representation 
matrices of SU{2) are, up to a normalization factor, orthogonal with respect to the Haar measure 
dfiH on SU{2): 

j dpi H (h)[7r f (h)] [^(h)] = ^r—6 fj 6 m r m 6 n , n (2.15) 

JSU(2) 47 + 1 

Moreover they serve as an orthonormal basis on the space L 2 (SU(2),dfiH) of functions f(h), h G 
SU(2), square integrable with respect to the Haar measure dfin on SU{2). Using the multilabels 

jj := (ji) • • • iJn), m '■= ( m ii ■ ■ ■ 5 m Jv) n '■= ( n i) • • • ? n 7v) j one can then define the so called Spin 
Network Functions (SNF) T j^-- for a graph 7 as 

T^ H (AU):= V / 2^TTK(M^I 7 ))] merie (2-16) 

e£E(~f) 

which serve then as an orthonormal basis on the Hilbert space 7i 7 = L 2 (A\ 1 , dfi^) by applying (12.151) 
for every copy of SU{2) belonging to the N edges (e\, . . . , ejv) in the edge set £'(7) of the graph 
7. Since every L 2 function / = f 7 (h ei , . . . , h eN ) o p 7 cylindrical with respect to a graph 7 can be 
expanded in terms of the spin network basis (I2.16P as 

j fa n 

with expansion coefficients zq A B \, it suffices to work with the basis states T j Aji (A\y) in H 7 = 
L 2 (A\~[, d/jLy). 

The inner product in Ti- can be extended to an inner product between SNFs TV?-,-., T' -,, 
belonging to different graphs 7, 7' by choosing a bigger graph 7" containing 7, 7' as proper subsets, 

16 Note that an su(2)-index can be shifted from upper to lower position (and vice versa) using the Cartan-Killing- 
metric Sij. However since 5ij = 5 y = ly is the unit matrix, the index position makes no difference. Hence we will 
adjust this position as necessary for notational convenience. 
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e.g. 7" = 7 U 7'0 Then one may extend the SNFs T^ Ait , T',-,,^,^, to depend on the additional edges 
according to the trivial representation, for example we set j e = in the extended version Ty,j AH of 
T j^j. if e G {E(y ) \ £'(7)}. In this way one arrives at the following inner product: 

{■^/mn'-^y'j'm'n') = ^7 7' f $m m' $n n' (2.18) 



2.2.3 Kinematical Hilbert Space Hkin of Loop Quantum Gravity 

The idea is then to construct the full kinematical Hilbert space Hkin '■= L 2 (A, dfio) of Loop Quantum 
Gravity, formulated for the complete uncountable set of all graphs (which can be embedded in E) as 
an inductive limit Hilbert space Hkin of directed sets of Hilbert spaces . . . C 7i 7 C Hy C Hy C . . . 
formulated over (infinite) directed sets of graphs . . . C 7 C 7' C 7" C . . ., where the direction goes 
to finer and finer graphs. Hkin can then be thought of as containing every possible such (infinite) 
directed set. Using (12.111) the measure dfio on Hkin, called the Ashtekar-Lewandowski measure, is 
constructed such that for the measure fio(f) of a function / e Hkin with / = / 7 op 7 one has: 

M/) : = / fdfJ'O = J (f-yop^)dfx = J / 7 ((p 7 )*d// ) = / f^dfh (2-19) 

that is cfyx 7 can be written as a push forward (p 7 )* dfi and we have the map p 7 : Hkin ~> Hy. 

Note that choosing the SNFs as a basis for Hkin, and demanding the validity of ( 12.181) for every 
pair 7 7^ 7', comes with some technical subtleties resulting from the possibility of splitting an edge. 
Consider a graph 7 containing only one edge e = £'(7). Now we may introduce an arbitrary number K 
of 2-valent vertices resulting in a formally different graph 7 7^ 7 with edge set E(j) = {ei, e 2 , . . . , ex} 
with e = ei o e 2 o . . . o e K . Therefore instead of having T^j mn {A) := [^-(^(A))] , we then find a 

new set of SNFs formulated with respect to 7: T-j^^A) = Ylk=i \ 7i jkQ l e k {A))~\ m , with quantum 
numbers jk = j, m k 7^ m,n k 7^ n. This process of splitting one edge can in principle be used to create 
arbitrarily many graphs 7 from one edge. However these different SNFs are (although formulated 
with respect to different graphs 7 7^ 7) not necessarily even orthonormal, and the question arises as to 
which subdivision of an edge one should initially take. Either one has restrict oneself to working with 
a fixed graph 7, or one has to impose additional conditions on how edges are split and/or composed. 
It turns out that imposing gauge invariance will solve this technical subtlety as a side-effect. 

The inductive limit construction of H k i n comes with the benefit that in actual computations one 
can restrict oneself to the finite dimensional graph level 7i 7 . As long as the structure of the directed 
set is respected (this is referred to as cylindrical consistency), the results obtained on the level of 
a single graph can then be extended to Hkin in an obvious way. In what follows we will restrict 
ourselves to the level of a single graph, that is the finite dimensional Hilbert space H^. 



2.2 .4 Gauge Invariance 

As can be seen in PQ, the Gauss constraints Cj as defined in ( 12. 6ft form an ideal of the constraint alge- 
bra. Therefore they can be solved separately from the other constraints (12. 4p . (12. 5p by constructing 
gauge invariant spin network basis states. 

One starts with a spin network basis state T a^J^A) as defined in (I2.16p . Using the possibility 
of redirecting and splitting edges in the edge set £(7) of the underlying graph 7 mentioned in 

17 Of course one has to be careful in cases where edges of 7, 7' intersect. Then one has to introduce additional 
vertices at these intersection points and to subdivide the intersecting edges, but this extension is straightforward. 
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the previous subsection, as well as the representation properties of the edge holonomies^j, we can 
introduce a new edge set £"(7), where each vertex v in the vertex set V(t) has only outgoing edges. 
Then we can rewrite T 



(A) as 



v /2jFT%,(/i(A))] -M 



(2.20) 



where we take the tensor product now over the set of edges E v {pj) outgoing at each vertex v G V(t), 
and the matrix M E ^ contracts all representation matrices {h(A))\ , fa- (h(A))\ of two edges 
ei,e2 if their composition gives an edge e of the original edge set £(7). (For example e\ o (e^) -1 — e, 
which implies [ir je (h e (A))] mn = £ r [tt,- (^(A))] mr [(tt^- (/i a (A))) -1 ] m .) Moreover, in such cases, 

Jei — Je 2 = je- in this way we have introduced trivial two-valent vertices vi trm ^ on each edge 
e G -£(7), which do not change the gauge behavior of the original function T j Aii (A) under a local 
gauge transformation g G Q, due to the transformation behavior (12.91) of the holonomy. 

Now by construction all edges in E(j) start at a vertex and end at a trivial vertex Ve , at 
which the gauge transformation glve™^) does not have any effect by construction. (Since at every 



trivial vertex v e the representation ^Kj e (g(v e nv ))] of the local gauge transformation g 
multiplied by its inverse due to (12.91) .) It follows that 



(triv) \ 



IS 



y/2fr+i[*Mv))] ■ M'M)] -M E(7) 



v t)6V(7) e6E„(7) 



(2.21) 



Now for four matrices (A, B), (C, D), each pair consisting of square matrices of the same dimension, 
we have (A • B) ® (C ■ D) = (A® C) • (B ® D). We can then rewrite as 



[*•*(*(«))] 



v «eV(7) ee-E„(7) 



M 



E(y) 



(2.22) 



It is obvious that we can obtain gauge invariant states T^jf(A) if we use the fact that the ten- 
sor product of two irreducible representations tt^,^ of 577(2) can be reduced to a direct sum 
®ji2=b'i-js| ^J' 12 °^ inducible representations 7Tj 12 whose matrix elements are constructed by using 
certain (anti-)symmetrizations of tensor product matrix elements [52l [56] . Therefore we can reduce 
the tensor product of representation matrix element functions at each vertex into a direct sum. 
Choosing a particular reduction order we can start with 



MO] 



ei,e 2 ,...,ejv„£-E„(7) 



K< 

3^2) = \h x -5z 2 \ 



MO] 



e3,...,ejv„6So(7) 



(2.23) 



where we denote the valence of a particular vertex v by N v . By successive reduction the whole tensor 
product of representation matrix elements can be decomposed into a direct sum 



ei,e2,...,ejv„G-B„(7) 



a 2 (jej 3e 2 ) a 3 ( a 2 jgg ) J (<W„-1 jVjy ) 



(2.24) 



*That is M/i e AJ] mn = Er h ft* )] mr h ft 2 )] r „ and [^(fe- 1 )]^ = [(7r 3 -(M) '] 
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The trivial representations (J = 0) contained in (I2.24p then provide gauge invariant states, that is 
we can introduce intertwiners / = {I v } ve v(-y) which project onto the trivial representations at each 
vertex v G V(^y) which are contained in the decomposition (12.241) : 



eG£„(7) 



M 



Eh) 



(2.25) 



To summarize: Taking certain (anti-)symmetrizations of spin network functions T j^^A), we can 
formulate at each vertex v G Vi^y) intertwiners projecting onto the trivial subspaces contained in 
the (anti-) symmetrizations. Each particular set / = {I v }v£V{-y) of projections then defines a gauge 
invariant state T^f(A). 

It can be shown that the T jf{A) provide an orthonormal basis of the gauge invariant Hilbert 
space Ti. = L 2 (A/G,dfi ), that is 



T q,p , T j f ) - 5y 7 5j,j 6pj 



(2.26) 



Note that the measure dfio is already gauge invariant, due to the translation invariance of the Haar 
measure dfin on SU(2). 

When gauge invariance is imposed the subtlety on possible subdivision of edges mentioned at the 
end of the previous section goes away: Consider again a graph 7 consisting only of one closed 
loop e, its beginning/end point being the only vertex of 7. It comes with a SNF T 1 ji{A) : = 
tij [7Tj(/i e (A))] . We now subdivide e into two segments e = t\ o e2, introducing a new graph 7 which 

has a trivial two-valent vertex vf . We can then formally introduce a set of SNFs T~^ n (A) : = 
[7r i (/i ei (A))] mini [7r i (/i e2 (A))] m2n2 . However gauge invariance enforces % f (A) := tr^ ^(h^A) ■ 

h e2 (A))~\ = Tyji(A). That is, spin network functions formulated over arbitrary splits of one edge are 
identified upon imposing gauge invariance. 



2.2.5 Spin States 

As it turns out, the action ( 12. 14ft of the right invariant vector fields on spin network functions 
TjjriiH can De re l a ted to the action of angular momentum operators Jg as used in ordinary quantum 



mechanics by [191 [50 1 [52]: 



-l) ]e+me Je I je 



-m e ;n e 



> 



(2.27) 



where | j e — m e ; n e ) denotes a usual angular momentum state | j e — m e ) with an additional 
quantum number n e which is not affected by the action of the angular momentum operator j£ . 

This opens the possibility of reformulating the theory using techniques from recoupling theory of 
angular momenta. In particular we may replace the tensor products at the vertices in (12.251) 



y/2h+l[ 



e€E v (j) 



eeE v (j) 



Yl I h mg- ; Tig" ) 



(2.28) 



We may express the gauge invariant intertwiner expressions in terms of so called recoupling scheme 
standard basis stated |a J — M — ; j n v ^,y Each of these states corresponds to a particular 



19 



See the appendix of the companion paper |52j for details and definitions. 
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trivial (J = 0) representation contained in the reduction (12.24)) . Fix a labeling of the N v edges 
contained in the set |ei, . . . , cn v } = E v (y). Each of the trivial representations contained in (12.241) is 
then uniquely labelled by the list of weights a 2 , a 3 , . . . , a at -i of the intermediate reduction steps in 
(E2I): 



I,, ■ 



?e£„(7) 



yj2je + i[iCj-(h(A))] = |a 2 (ii,j2) a 3 (a2,j 3 ) . . . a,N v -i(aN v -2, jjv„-i) J(a w „-i, ji\r„) = M = ; j v ti v ) 



such that 



|o J=0 M = ; j„ n„). 



(2.29) 



T — 



\a J = M = ; j v n v ) 



(v) 



■M, 



E( 7 ) 



(2.30) 



See figure [H for an illustration. 



3 The Volume Operator 

3.1 Definition of the Volume Operator 

The operator corresponding to the volume V(R) of a spatial region R C £ 



V(fl) = / d 3 x y/detq(x) = / d 3 a; 



acting on gauge invariant] 20 ! spin network states, as derived in [40, 52J is defined as 



(3.1) 



d pjdet(q(p)) 1 T 



d 3 p V» 7 T* 



where 



^,7 



i<j<k<n 
(Jij) 2 , (Jjk) 



(3.2) 



(3.3) 

(3.4) 
(3.5) 



and (J/j) 2 = ELi( J / + Jk jf - The construction of (E3D, (E3J), (ESI) is explained in detail in [52]. 
Here £p denotes the Planck length, and Z is a real constant whose specific value was found in [19] as 
^ _ p33hi__ p or our p Ur p 0ses we se t Z to 1 or i, since the Immirzi parameter (3 can be freely chosen 
and Z will only give an overall (at every vertex identical) constant scaling. The choice of Z = 1 vs. 



20 We restrict ourselves to the gauge invariant sector in the following. Almost everything could easily be formulated 
in terms of gauge variant spin network states. 
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i has the effect of rotating the spectrum in the complex plane from the imaginary axis (Z = 1) onto 
the real axis (Z — i), which clearly makes no difference. 

The sum in (13. 3 p has to be taken over all vertices v G V(7) (each vertex having N v out- 
going edges) of the underlying graph 7. At each vertex v one has to sum over all possible or- 
dered triples (ei,ej,ex) = (UK), I < J < K < N v of outgoing edges at v. Here e(IJK) = 
sgn( det (ej(v ), e j(v), ck(v ))) denotes the sign of the determinant of the tangents of the three edges 
e/, ej, ex intersecting at v. 

As is obvious from (13.31) . the volume operator acts nontriviallv only at the vertices of the graph 7. 
Moreover its action at distinct vertices is completely unrelated 21 !. We can therefore restrict ourselves 
to analyze (13. Ah at single gauge invariant vertices, since every graph 7 and hence every T jf can be 
constructed from those as can be seen from (12.301) . In what follows we will therefore consider only a 
single vertex v and hence the single vertex operator V := V va . Moreover we will denote the valence 
of the vertex by N, dropping the label v. 



3.2 Computational Task 
3.2.1 Matrix Formulation 

The remaining task is then to consider 



V 



z . < IJK ) 



Qijk 



KJ<K<N 



ZQ 



Q ] Q 



(3-6) 



where Q is a totally antisymmetric matrix with purely real elements. Its eigenvalues Aq are purely 

imaginary and come in pairs Aq = ±i • A, A > 0, A 6 R@ We choose Z — 1. Equation (I3.6P then 

1 

tells us that the volume operator V has the same eigenstates as Q and its eigenvalues A^ = |Aq| 2 . 
Thus we are left with the task of calculating the spectra of totally antisymmetric real matrices of 
the form: 

Q:= <UK)qijK (3.7) 

KJ<K<N 



3.2.2 Computational Realization of the Recoupling Scheme Basis 

Here we briefly comment on the computation of the recoupling basis. In some cases the number of 
gauge invariant states contained in all possible recouplings of a set of spins can be surprisingly small. 
We are not aware of a general exact analytic treatment of this problem, only the case where all spins 
are taken to be j = \ [59J. 

A standard recoupling scheme is the successive recoupling of all iV angular momenta ji, . . . , Jn 
at the vertex v, as follows. Fix a labeling of the edges. Then recouple ji,j2 to a resulting angular 
momentum 02(31 32)- We have, due to the Clebsch-Gordan-Theorem, |ji — j 2 1 < 02(31 32) < ji +32- 
Now couple 0.2, J3 to a resulting angular momentum 03(02 J3) where again \a<z — j%\ < 03(02 33) < 
a 2 + jz- This continues until we arrive at the last recoupling step where J = a N (a N ^x j N ) and thus 

21 Though of course the vertices are related through the linking structure of the underlying graph 7. 
22 Q^Q is a totally symmetric real matrix, whose eigenvalues AgtQ are purely real and come in pairs AgtQ = |Aq| . 
This provides an alternate approach to computing the volume eigenvalues, which is discussed in [52) . 
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|ajv— i — 3n\ < J < a N~i + Jjv- If we work with gauge invariant recoupling schemes we have J = 

and therefore we must have ]n = a^-i- 

The challenge is now to find, for a given spin configuration ji, . . . all possible gauge invari- 



ant states, that is all recoupling sequences as illustrated below for which a k 



(min) 



K-i - jk\ < 



ak{ak-i 3k) < (ofc-i +j k ) 



a 



(max) 



(with ai = ji) and additionally a n _i 



]n is fulfilled [52] 



J 



a 



(max) 



(max) 

(max) 



(max) 
0>2 



(min) 



II, 



(max) 







"min) 



(min) 



min) 



Jmin) 
A 6 



a 2 



0-3 04 a 5 a 6 



a 



(max) 
N-3 



a 



(max) 
N-2 



a 



(max) 
N-l 



(min) 
l N-3 



(I 



(min) 
N-2 
(min) 
N-T~ 







a 



(max) 
N 



OTV-3 O-N-2 a N-l 



J = 



Figure 1: Recoupling Scheme: A schematic recoupling of N spins, labeled by the inter- 
mediate steps 02, 03, • • •• The final recoupling step must go to total angular momentum 
zero. c.f. section 13.2.21 



3.2.3 Gauge Invariance at the Operator Level 

At the operator level gauge invariance 

3 N 



k=l K=l 



implies that 



N-l 



Jn = -J2 J l Vfc=l,2,3 



(3.* 



(3.9) 



L=l 



This can be used in ( 13. 6j) in order to eliminate all <jun containing N as an edge label [52]. Hence 
the volume operator can be rewritten as 



V 



QlJK 

KJ<K<N 



(3.10) 
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where 



a(IJK) := e(IJK) - e(IJN) + e(IKN) - e(JKN) 



(3.11) 



This is the final expression for the volume operator acting on gauge invariant vertex recoupling states. 
Expression (13.101) generalizes the observation for the special case of a 4-valent vertex as stated in 
P] . section V.D, and 0]. 



3.3 Implementation 

Obviously the analysis of (13.10!) contains two parts: a recoupling computation of the matrices qux 
and a combinatorial analysis of the possible sign configurations a(IJK). 



3.3.1 The quK constituent matrices 

For a given spin configuration jx, . . . , at an N— valent vertex v the trivial representations contained 
in the tensor product decomposition (I2.24p . that is the gauge invariant recoupling standard basis 
|aJ = 0M = 0; n ) as given in (I2.29p . can be computed according to section [3.2.21 A simplified 
general expression for the matrix elements of qux was derived in [51] and will not be reproduced 
here. An even more detailed discussion^! of this expression necessary for the present analysis is 
presented in [52] . 

As we have stated before, an arbitrary but fixed labeling of the edges has to be chosen in order to 
perform the reduction (I2.24p . Since a relabeling of the edges can be implemented as a unitary basis 
transformation among the recoupling states (12.291) . the spectrum of the volume operator is invariant 
under edge permutations. (See [52] for a more detailed discussion.) We can use this property to 
drastically decrease the number of spin configurations to consider in the numerical analysis, because 
we can always choose a particular edge labeling ei, . . . , in which the spins are sorted 

jl < h < • • • < 3N = jmax (3.12) 

Each of these spin configurations has degeneracy D(j) as follows. Assume we have among the N spins 
j = (ji, . . . , ]n) V subsets of identical spins, each such subset containing Nk (k = 1 . . .p) elements. 
Such a configuration then stands for 

D(J) = W .N-A-...-N r > (3 - 13) 

spin configurations. Moreover we have the conditions J2n=i 3n > 3n and ^ n=1 j n G N which are 
necessary in order to obtain a trivial representation in the reduction (I2.24p . 



3.3.2 Sign Configurations 

In [57] it was shown that for vertices of valence greater than 4, which are formed by self inter- 
secting^] smooth edges, there exist continuous moduli parameters labeling the thus formed knot in 
a diffeomorphism invariant way. However there is also a discrete set of information regarding the 
linear dependence of edge tangents contained in the tangent space at the intersection point: Three 
tangents can be coplanar or linearly independent. Moreover, as we are going to show now, there 

23 There exist certain special edge combinations which have to be worked out carefully. 
24 These are essentially embeddings of S 1 — > S. 
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is diffeomorphism invariant information on the relative orientation (right handed / left handed) of 
linearly independent triples of tangents: If we have five or more edges at a vertex, then we can look 
at the orientation of each triple of tangents at that vertex. We find that (up to an overall orientation 
flip) only special combinations of tangent orientations can be embedded in E. 

Let us denote by e the set of signs e(IJK) = 0, ±1, one for each ordered edge triple (ej, ej, ex) 
with I < J < K < N. We have (£) such signs describing a particular embedding of an iV-valent 
vertex v. According to this definition we denote by a the set of values a (UK) = 0, ±1, ±2, ±3, ±4 
computed from a particular set e via (13.111) for each of the { N ^ X ) ordered edge triples (ej, ej, e#) 
with I < J < K < N . Note that the map S : e i— > a is non-injective, in that there can be 
more than one ^-configuration giving the same c?-configuration. We will denote this multiplicity by 
Xs '■— | {sign configurations e for which : S(e) — a}\. 

Because of the complicated structure of the recoupling computation, to our knowledge no one has 
yet attempted a detailed combinatorial analysis of the edge tangent information for vertices of valence 
greater than 4 and its effect on the volume spectrum. In particular it is not even obvious which triple 
sign configurations e can be realized at a vertex embedded in a three dimensional hypersurface S 
and what a combinations result from that in the gauge invariant sector. 

In the present analysis we have computed the e configurations by a Monte-Carlo process, as 
follows. Consider an iV-valent vertex v. Now randomly sprinkle iV (labeled) points on the surface 
of a unit sphere and compute the sign e(IJK) for each of the ( 3 ) triples (UK). This gives a e- 
configuration, which is stored in a (hash) table, along with the number of times that e-configuration 
has arisen so far in the sprinkling. At the beginning of the sprinkling process there will always be new 
sign configurations obtained which are successively added to the table. However, after some time, all 
the obtained sign configurations will already be in the table. This can be taken as an indication that 
all possible sign configurations have been obtained, so the computation is terminated. (In particular 
we stop sprinkling when the smallest number of 'hits' associated to a e-configuration is greater than 
one.) 

A disadvantage of the sprinkling process is that one cannot be sure to have obtained every possible 
sign configuration. Furthermore, the probability of sprinkling three points onto a single plane is zero, 
so the sprinkling yields only non-coplanar sign configurations. This is the reason why in the present 
analysis we consider only vertices where each triple of edge tangents is linearly independent. We will 
continue to call vertices with such edge configurations non-coplanar vertices. 

Note that for non-coplanar vertices we have e(IJK) = ±1 for all triples (UK) and therefore each 
a(IJK) may assume, due to (13.111) . the 5 values a(IJK) = 0, ±2, ±4. Thus there are in principle 

5^ 3 i different a-configurations possible. The number of (^-configurations given here is the number 
of numerically distinct a-configurations which occur when imposing gauge invariance on each of the 
realized e-configurations. The number of a — configurations (i.e. those for which each of the 
( N ^ 1 ) cr(IJK) are zero) is given by the number of occurring e-configurations which lead to a a = 
configuration. The sprinkling process gives the following results: 

Due to recent progress it now seems possible to overcome the mentioned disadvantages of the 
sprinkling process, as it is feasible to compute the possible e-sign configurations directly, without 
using Monte-Carlo methods. The idea behind this computation is given in |52j and will be worked 
out in more detail in a forthcoming paper [53]. Up to now we have verified that the number of 
realized e*-sign configurations of a completely non-coplanar 5 & 6-valent vertex is indeed given by 
exactly 384 and 23,808 respectively. 
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Table 1: Sign factor combinatorics for 4-7-valent non-coplanar vertices 



We would like to mention here another interesting implication of the sign combinatorics given in 
the last column of the table above. Equation (I3.10P can be interpreted as a kind of 'self-regulating' 
property of the volume operator when applied to certain 'pathological' edge configuration^! [52"] . For 
vertices of valence greater than 3 there exist vertex geometries containing only non-zero edge triples, 
however all of them occur with identical sign, e.g. e(IJK) = 1 V / < J < K. Upon demanding 
gauge invariance such vertex geometries will have zero volume due to (13.101) and (13. lip , a statement 
which is independent of the spin configuration at a vertex@ 

4 Spectral Analysis: Numerical Studies 

With the preparation of the previous section we now have all prerequisites for an analysis of the 
spectral properties of the volume operator. While our computational implementation is applicable 
to arbitrary vertex configurations, it is limited in practice, of course, by finite computer resources. 

4.1 Computer Resources 

We performed most computations on the Peyote high performance dual 64-bit 3.4 GHz Xeon pro- 
cessor cluster, at the Max Planck Institute for Gravitational Physics in Potsdam, Germany, using 
perhaps 20,000 CPU-hours in total. Some computations were also performed on a 2.2 GHz dual 
Opteron machine at Imperial College. The eigenvalues were stored in 188GB worth of binary data 
files. 

4.1.1 Code and Speed Issues 

Our code is implemented within the Cactus framework (described below), and employs the LAPACK 
library to compute eigenvalues using singular value decomposition, the details of which are described 
in [52]. The set of cr-configurations for which we compute eigenvalues is distributed among the 
processors. 

25 In [25 an upper bound for the expectation value of the operator corresponding to the inverse of the classical scale 
factor in cosmological models is given and found to be finite. However, this bound grows with the number of edges at 
a vertex, and an explicit vertex configuration is given there for which this indeed happens. However, this configuration 
gives zero volume upon imposing gauge invariance, so it seems feasible that the problem with the explicit valence 
dependence of the upper bound may be cured. 

26 0f course the spin configuration has to allow for gauge invariance. 
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Because of the smallest non-zero eigenvalues which decrease with j max , our code is in fact bounded 
by the 64 bit numerical precision for the 5- vertex (and 6- vertex to a tiny extent.) We can carry the 
machine computations beyond j max = y, but the smallest non-zero eigenvalues will descend deeper 
into the numerical noise. See section 14.21 At higher valences our computations were bound in effect 
by disk space, because of the enormous number of <r-configurations at large valence. In order to 
test our code we also performed computations on the 4-vertex. These computations were bounded 
simply because we use a single byte to store the edge spins. We could go much further than the 
computations in [ST], if desired. 



4.1.2 The Cactus Framework 

Integration of our code into the Cactus high performance computing framework [61] affords a number 
of benefits, which is discussed more fully in [52]. In summary, the Cactus framework provides 
automatic parallelism, facilitates linking to LAPACK and Fortran code, and allows a great degree 
of portability. We hope that the true value of implementing our code within the Cactus framework, 
however, will be that it may now be readily used by other groups, as a prototype for performing 
numerical Loop Quantum Gravity computations, for the full theory. Instructions on how to download 
the code used for the computations described in the paper will appear shortly. The intent is to begin 
developing a community code base for such computations, taking advantage of the modularity and 
collaborative nature of Cactus to facilitate further numerical investigations. 



4.1.3 Actual Computation 

We present in this section some results of a numerical study for 5-7-valent vertices, where we have 
computed the matrix representations of the volume operator for all sets of edge spins j± < . . . < 
3n = jmax (A = 5, 6, 7) and for all a-configurations according to table [TJ As we want to treat all 
^-configurations for a given valence, we can take j max up to a certain cutoff which is strongly limited 
in practice by the number of cx-configurations. 

The eigenvalues of the computed matrices are obtained using singular value decomposition, which 
is applicable due to the modulus in the definition (13.101) of V. The eigenvalues are finally collected 
into histograms. The range of eigenvalues is divided into a certain number Nuns of intervals. Then 
for each such interval, which we refer to as a histogram bin, the number of eigenvalues falling into it 
(occupation number) is counted. Here each such eigenvalue is regarded as coming with a redundancy 
D(j) ■ Xs where D(j) is given in (13. 13j) and as described above Xff denotes the number of e-sign 
configurations giving the particular ^-configuration. The occupation number of the histogram bins 
then serves as a notion of spectral density. We can discuss the properties of this spectral density, 
when superposing the eigenvalues coming from different parameter ranges. We omit the eigenvalues 
Xy = from all histograms, as they are much more numerous than any of the displayed occupancy 
numbers. 



Table [5] illustrates different ways to obtain histograms. For example, case 3 shows that for a 
fixed vertex valence N we can collect all eigenvalues from spin configurations j taking all values up 
to the spin- value cutoff j max , for fixed a, in order to see the effect of different vertex embeddings on 



the spectrum. We may also, as in | case 4 1 , combine all fixed a histograms for a given vertex, but 
restrict to a fixed j max = j^, to illustrate the dependence of the spectrum on j max . Finally, as shown 
in 



case 6 , we may take all eigenvalues we have computed into one histogram. In each histogram 
we plot the occupation number iV evals of the bins versus the eigenvalue X v . Note that the displayed 
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case 


valence N 


tr-conng 


Jmax — JN 


j-config 


1 


fixed 


fixed 


fixed 


fixed 


2 


fixed 


fixed 


fixed 


all 


3 


fixed 


fixed 


all 


all 


4 


fixed 


all 


fixed 


all 


5 


fixed 


all 


all 


all 


6 


all 


all 


all 


all 



Table 2: Some possible choices of parameter ranges for histogram construction: 'fixed' implies keeping 
the according parameter fixed, 'all' implies placing all eigenvalues in the histogram, regardless of the 
value of this parameter. 



eigenvalues \ v have to be multiplied by y/\Z\£p, where due to our conventions (13.41) . (13.10p . Z is a 
constant and £p is the Planck length cubed. 



4.2 Results 

Our numerical data reveals that for vertex valence N > 4 the sign factors contained in the defini- 
tion (13.41) of the Ashtekar-Lewandowski volume operator open a fascinating interplay between spin 
quantum numbers and gauge invariance on the one hand, and the graph label on the other hand, 
which manifests itself in different sign- dependent spectral densities of this operator. This has not 
been analyzed in detail before due to the already very complicated recoupling computation. 

In general we find that the actual form of an eigenvalue distribution heavily depends on the 
sign factor combination a which is induced by the vertex geometry ( case 3 in table [2]), whereas 
varying j max while taking all a-configurations ( case 4 in table [2D only seems to rescale the eigen- 
value distribution, since a larger j max gives more eigenvalues. This is illustrated in figures [5] and [3] 
respectively, for the gauge invariant 5-valent vertex, which we have computed up to j max = y. As 
shown in tabled] we have 171 a-configurations for the 5- vertex. Due to the modulus in the definition 
(13.101) of V two sign configurations a, a' give the same spectrum if they differ by an overall sign only, 
that is a(IJK) = —cr'(IJK) W, J, K. Hence (excluding a = ) we effectively have 85 distinct o 
configurations for the 5- vertex. 

Remarkably the smallest non-zero eigenvalues in a dataset can behave in three different ways 
as the largest spin j max at the vertex is taken to higher values, depending on the (^-configuration 
at the vertex. They can increase (as was already known for the 4- vertex), but they can also stay 
constant or decrease. For the latter behavior we find sequences of smallest eigenvalues which tend 
to decay towards zero exponentially with j max , for the 5- and 6- vertex. For the 7- vertex this is 
difficult to decide, as we do not have enough data points, due to limited computational resources. 
We however find all three behaviors of the smallest non-zero eigenvalues for the 5, 6, and 7- vertex. 
This is illustrated in figure HI 

In the case of the histograms for the 5- vertex this results in an accumulation of eigenvalues close 
to zero, as can be seen in the inset plot of figure El When going to larger eigenvalues the spectral 
density decreases to a minimum and then rises again. This second rising edge can be fitted with an 
exponential, which one would suppose by demanding a quasi-continuous volume spectrum at large 
spins. The gauge invariant 6-vertex, shown in figure [51 also possesses the rising edge, but there is 
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1 .4e+07 




450 500 



Figure 2: Gauge invariant 5- vertex, case 3 of table [2) Histogram for each of the 85 a- configurations. 
For each fixed a all eigenvalues resulting from spin configurations j up to j max = y are collected. 

Note that we have treated all ^-configurations equally here (xff = 1) in order to make a direct 
comparison possible. All but ten of the histograms have iVbi ns = 2048, the others N^, ins = 512. These 
numbers are chosen to generate the sharpest curves possible. 



9e+09 




Figure 3: Gauge invariant 5- vertex, case 4 of table [2) Histograms up to j max = y (top curve, below 
it are histograms for smaller j max )- N hins = 2048. For j max = y there are 4.6 x 10 12 eigenvalues in 
all. Inset: Portion of histogram for Xy < 9. 
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Figure 4: Smallest non-zero eigenvalue A^™ for each of the 85 a-configurations (arbitrarily labelled) 
as a function of j max , for the gauge invariant 5-vertex. Note that multiple <r-configurations lead to 
identical A^™ m ^ (j max ) . This is presented in detail in the companion paper [52] . 



1.8e+10 




20 40 60 80 100 120 

Figure 5: Gauge invariant 6- vertex, case 4 of table [2j Histograms up to j max = -y (top curve, below 
it are histograms for smaller j max ). N^, ins = 2048. There are 7.4 x 10 12 eigenvalues in all. 
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no minimum at small eigenvalues in the spectral density, the eigenvalues seem to be rather equally 
distributed close to zeroo For the gauge invariant 7-vertex we obtain a rising edge similar to the 
6- vertex, without a lip. However as can be seen in [22], it is not obvious if the spectral density rises 
exponentially. This might be caused by the fact that, due to the huge number of cx-configurations 
at the 7-vertex, we have calculated the spectra only up to j max = |, so the spectrum seems to be 
largely dominated by the proliferation of cx-configurations. Due to the necessary cutoff for j max the 
rising edge in all histograms then merges into a decreasing tail (which can be considered entirely a 
cutoff effect), however the rising edge can be extended towards larger eigenvalues as the cutoff for 
jmax is increased. 

The behavior of the largest eigenvalues A^ nax ' 1 also depends on the ex-configuration, as shown in 
figure [6] for the 5- vertex. 




Figure 6: Largest eigenvalues vs. j max for each of the 85 ex-configurations of the gauge invariant 
5- vertex. Note that many ex-configurations lead to identical X^ ax \j max ). 



Figures H] and El suggest the presence of a non-trivial symmetry of the volume operator. In partic- 
ular, they show groups of ex-configurations possessing identical smallest or largest eigenvalues. These 
ex-configurations can be grouped into equivalence classes, which is done explicitly in [52J. Unlike the 
ex — > — ex symmetry, it is not obvious why such ex-configurations yield identical volume eigenvalues. 
One might presume that this is caused by the fact that a permutation of the edge labels at a vertex 
acts as a unitary transformation (see section 13.3.11) on the volume operator and in fact one has to 



This absence of a non-zero minimum in the spectral density at the 6- vertex is not necessarily maintained when one 
restricts the histogram to a single a-configuration, as one can see in figure 25 of [55]. Thus it may be the 'averaging' 
over (7-configurations that is removing the 'lip' at zero, which may reemerge only at larger j max (c.f. section 6.3.1 of 

152] )■ 
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look at the orbits of the a-configurations under the permutation group [53] together with a possible 
presence of groups of identical edge spins contained in spin configurations contributing the small- 
est / largest eigenvalues. 

The cubic 6-vertex is used often, because from it one can easily generate a graph which forms 
a cubic lattice. Although it possesses coplanar edge triples, we have also analyzed it here, and find 
that it belongs to the class of a- configurations whose smallest non-zero eigenvalue rises as j max is 
increased. 

As one can see from equation (13.101) . the number of a(IJK) rises with the vertex valence N as 
( Hence to compute the complete matrix representation of V we have to add the according 



number of constituent matrices qjjK- Figure [7] shows an overall histogram ( case 6 in table [2]), 
where all eigenvalues we have computed are depicted for valences N= 4-7. Note that the effect of 
the rising number of edge triples is obvious in figure [7J for example a 7-valent vertex can contribute 
relatively large volume eigenvalues from relatively small spins. This seems to be a general feature, 
that very high valent vertices can produce very large^l volume eigenvalues, contributed by very small 
spins. Such spectra are expected to be dictated by the edge combinatorics at the vertex, rather than 
the recoupling of spins. That is, this would give rise to 'macroscopic quantum objects', as is the 
case when one observes quantum mechanical behavior not only for single carbon atoms but also for 
C 60 -molecules. 



1e+12 



1e+10 



1e+08 



4- vertex 
4, 5- vertices 
4, 5, 6- vertices 
4,5,6, 7-vcrticcs " 



10000 
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10000 



Figure 7: Histogram of entire dataset, in log-log scale. iVbins = 4096. Each successive curve adds 

™ for the 4- vertex, y 



eigenvalues for the next larger valence. Here the cutoffs for j max are: hp for the 4-vertex, -~ for the 



5- vertex, y for the 6-vertex, and | for the 7- vertex. 

We find the following characteristic spectral parameters for the different valences: The scaling 



i Maybe even macroscopic if we go to very high valences. 
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of the overall (all cx-configurations) smallest eigenvalue sequence A^ 1 (jmax), the overall (all cr- 

configurations) largest eigenvalue sequence X^ ax \j m ax), the fitted number of eigenvalues N eva i s (Xy) 
at the rising edge of the histogram for fixed j max , and the total number of eigenvalues (including 
redundancies) A^'°2 (j max ) contained in a histogram. 



valence N 


A (min), s 


Umax),, s 
VJmax; 


^Vcvals(A^) 


N {tot) (i ) 

iv cvals Umax J 


5 


~ 30 • e- 79 A— 


° Umax,/ 


2.9 • 10 6 • e 43 A o 


~9.3-10 3 Umax ) 5 - 9 


6 


~ 40 • ce" 243 A— 


*~ 3.5 • (jmax) 1 * 1 


0.19 • 10 9 • ce - 55 A v 


-2.9-10 6 ( Jmax ) 7 - 5 


6 (cubic ) 




~6.5- Umax ) 145 






7 






- 0.37 • 10 9 • ce - 79 A v (?) 





Table 3: Fitted parameters of the eigenvalues and spectral densities for various vertices. 



Here we have chosen the fit functions by inspection of the data rather than from an analytic model. 
For example, we can plot the smallest eigenvalue sequence A^ mn ' ) ( j max ) in a log-log plot or a log-plot 

and then decide whether a polynomial dependence A^ 1 (jmax) ~ cl • (jmax) 6 or an exponential depen- 
dence A^ nm ' ) (j max ) ~ a ■ e 6 "3' ma * gives a better fit to the data. Note however that for the rising edge of 
the histograms one could expect an exponentially rising number of eigenvalues from the numerical 
analysis of the area operator [32] . The set of data points to be used for the fitting process is chosen 
manually by inspection of our data (see [52] for details). This should be kept in mind when looking at 
the numerical values of the fitted parameters given in table El Our point here is to demonstrate the 
possible behavior of the spectrum. For the 7-valent vertex the fit itself is questionable, as due to the 
huge number of cx-configurations our computation can only be extended to a rather small spin cutoff 
of |. Note that the growth of A^ nax ' ) (j max ) as ~ (jmax)^ can be understood from the actual form of 
the volume operator (13.101) . as it contains sums of g/j^-matrices, each coming from a polynomial of 
degree 3 in the spins, and finally the square root is taken. 

In [3B] a formula for the spectral density of the volume operator at the gauge invariant 4-valent 
vertex was derived in the limit of large spins, treating ji, . . . , j 4 as continuous classical variables. As 
can be seen from equation (15. ip below, the recoupling process of four angular momenta ji, . . . , j'4 to a 
resulting total angular momentum J = implies that ji, j 2 , j'3 have to be recoupled to a total angular 
momentum CI3 = j'4 := j. Now in [IE] two different eigenvalue densities for the matrix expression 
( 15. 3p of the volume operator are derived: Pj(Xy) and p(A^). The former gives an expression for the 
eigenvalue density Pj(^v)\jiJ2,j3 a ^ fixed j max = j 4 := j and fixed spins ji, j 2 , J3, which is what we 
would need here for comparison with our spectral density. We would have to integrate Pj(^v)\j 1 ,j 2 ,j3 
over all ji, j 2 , J3 in order to get a similar situation as we have considered in our histograms, where we 
collect for fixed j := j max all eigenvalues from spin configurations ji < j 2 < J3 < jmax- Unfortunately 
Pj(^y)\ji,j 2 ,j 3 integrated over ji,j 2 ,j3 cannot be accessed explicitly, as the integral over ji , J2 ? J3 
involves combinations of transcendental functions of the spins. What is accessible is a so called 
'integrated eigenvalue density' p(^v)h,j2,hi which is defined as p(Ay)|j 1 j 2 j 3 := / djPj(^v)\h,h,h- 
This is equivalent to considering for a fixed set ji,j 2 ,js all possible values of 03 = j 4 (see section 
I3.2.2p at once. This can be viewed as considering the spectral density of a direct sum of all matrices 
Q of (13. 6p for fixed ji,j 2 ,j3 over all j 4 = j, which can be reached by a recoupling of ji, j 2 , j s . While 
this is an interesting result, it is difficult to extract a notion of spectral density out of [IE] as we 
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have computed it throughout our numerical studies. However it might be worthwhile to try applying 
further approximation techniques in the integrations such that properties of the total eigenvalue 
distribution (integrated over all spins) can be extracted. 

Note on Numerical Errors 

As can be seen from table [3] there is a serious caveat when looking at the behavior of the smallest 
non-zero eigenvalues. As we go to larger j max these eigenvalues decay exponentially to zero. However 
as we increase j max the matrix entries in the matrices which are finally diagonalized become larger 
and larger. Hence when analyzing the smallest non-zero eigenvalues, we compute tiny numbers from 
relatively large numbers. (In fact it is this effect which limits how far we can carry the computations 
for the 5- vertex, as mentioned in section 14.1.11 ) 

In order to consider the classification of (^-configurations by the behavior of their smallest eigen- 
values as a true result of our analysis it is therefore mandatory to keep track of the numerical errors 
in our computations. This problem turns out to be technically very puzzling. We were able to gain 
a good distinction between numerical errors and true results and have verified the behaviors of the 
smallest eigenvalues above as true properties of the volume spectrum. A detailed discussion on this 
issue can be found in the companion paper [52] . 



5 Analytical Results on the Gauge Invariant 4- Vertex 



In this section we will present analytical results on the spectrum of the volume operator at the 
gauge invariant 4- vertex which will extend some of the earlier analysis in [51]. In particular we will 
give a short proof that the spectrum of the volume operator at a given gauge invariant 4-vertex is 
simple, that is all its eigenvalues (except zero) come in pairs, and there are no further degeneracies 
or accumulation points in the spectrum. Furthermore we will give an explicit expression for the 
eigenstates of the volume operator in terms of polynomials of its matrix elements and its eigenvalues. 
Thirdly we will state an explicit result on the scaling of the smallest non-zero eigenvalue of the volume 
operator with respect to the largest spin at the vertex. This will support the numerical results of 
[oTj . Note that the specific matrix realization of the volume operator plays a crucial role here: It can 
be written as an antisymmetric purely imaginary matrix having non-zero entries only on the first 
off-diagonal. This makes it possible to apply certain techniques from orthogonal polynomials and 
Jacobi-matrices [54J . This progress of course contains the drawback that in the presented results the 
structure of the matrix realization is heavily used. So the generalization to higher valent vertices will 
be technically very involved, if it is feasible at all. 

5.1 Setup 

The general case of the gauge invariant 4-vertex deals with 4 edges ei, . . . outgoing from the 
vertex v. Each edge (k = 1 . . .4) carries a 2jk + 1 dimensional representation of SU(2). Due to 
gauge invariance the standard recoupling basis states (12.291) are given in this 



where the intermediate recoupling 03(02^3) has to equal j'4 due to the Clebsch Gordan Theorem 
and as a result the intermediate recoupling 02 (ji 32) is the only degree of freedom for fixed spins 




(5.1) 
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ji, . . . , j'4. Therefore expression (I3.10p for the matrix representation of the volume operator simplifies 
dramatically to give 

V=yf\Z- a(123) q 123 \ = {J\ Z|2[a(123)] 2 • QtQ| (5.2) 

where Q := gi 23 and <x(123) = e(123)-e(124) + e(134)-e(234) = 0, ±2, ±4 gives a constant numerical 
prefactor depending on the relative orientations of the 4 edges only. In the following we will assume 
that <t(123) 7^ but leave its numerical value unspecified. Here we set Z — i. Using the explicit 
formula for the matrix element of the volume operator [5TJ [52] we obtain for the matrix Qd (D 
denotes the dimension) a matrix with the following shape 



/ 



Q 



D 



- qi 

qi 

q 2 

q 3 









-q 2 

-q 3 

-q 4 

<?4 -q 5 

q 5 













V 





























where the matrix elements are given by 

qk ~ v /(2a 2 +l)(2a 2 -l) 

X (j3 + h + «2 + 1) (-J3 + J4 + a 2 )(h ~ Ji + a 2 )(j3 + k -02 + 1) 
and the intermediate recoupling step a 2 (j\j 2 ) is given as a 2 = a^ 1 ^ + k, that is 
























































































■ 

-qD-5 
















QD-5 


-QD-4 













<?d-4 





-QD-3 













QD-3 





-qD 


-2 










QD-2 





-to- 


-l 
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l +32 + a 2 )(ji 


-h + a 2 )(ji +32 ~ 


- a 2 + 1) 


X 



max [|j 2 - |j 4 - 
and D = dimQ D = a^ 1 ~ o!^ 1 ' 11 ^ + 1. 



{rain) 



^ ^ (max) . r . . . . i 

<a 2 <a y 2 = mm [j 1 + j 2 , j 3 + j 4 J 



(5.3) 



(5.4) 



(5.5) 
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5.2 Non-Degeneracy of the Spectrum 

By inspection an eigenvector \I/ A = (\I> A , . . . , ^ D ) T for the eigenvalue A of (15. 3p with (Qd — Al)\I/ A = 
will fulfill: 



(I) -A^ A -gi^ A = 

(II) q k -^ti - ~ Qk*Ui = 
(III) gjv-i*k-i - A^ A = 



k = 1 








1 < fc < D 




k = D 





Now suppose A ^ 0@ Setting \1/ A = leads together with (I) to \& A = which then in turn 
implies together with (II) that \1/ A = V7c. Therefore in order to get a non-vanishing eigenvector \E fA , 

! 

A 7^ we must have \1/ A 7^ 0. 

Now suppose there exists an eigenvalue A 7^ of Qd that has multiplicity > 1 : Let \1/ A , \I /A be two 
eigenvectors for A. Then we can find two constants a, a 7^ 0, a, a G C such that we can construct 
an eigenvector $ A = a^/ x + a^/ x with the property $ A = 0. As we have seen this implies $ A = Vfc, 
and thus = a\l/ A + <3\E' A , that is \1/ A , \I/ A are linearly dependent and therefore there exists only one 
eigenvector \1/ A for each eigenvalue A 7^ of Qd- We have thus proved the following 

Theorem 5.1 Spectral Simplicity of Qd 

(i) The spectrum spec(Qo) ^ is simple (consists of D real numbers). 

(li) If^ x = (^ A , . . . , ^ D ) T is an eigenvector ofQ D (Qd^ x = A*\ A 7^ 0) then ^ A 7^ 0, ^ D 7^ 0. 



5.3 Eigenstates in Terms of Eigenvalues 



Starting from (I) in (15.61) by setting \I/ A = x = const one can explicitly construct the components 
\I/ A . One getqfj (using the integer number L > 1 with L < — for even D and L < for odd D): 



-xA 



2L 



2L-1 

n ?fc 

fc=i 



L-l 



M=l 



2L+1 - 2L 

n 9* 



A 2L + £ 6g L+1) ■ A 2 ( L " M > 



M=l 



fc=l 



where we have defined 



■,(*o) .. 



k M -2 fei-l-2 fc 3 -2fc 2 -2 

E E - E -EE Qk M QkM-i a .k 2 a .ki 

k M =2M-l fe Af _ 1 =2(M-l)-l fcj=2(Af-0-l fc 2 =3fc 1= l 



(5.7) 



(5.* 



29 This case has been already discussed in [51] , 

30 Note that (|5.7[) can be compared to results of [JS] which where communicated to the authors during the preparation 
of this paper. 
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Note that if M = 1 in this sum one has to take the upper bound k — 2 of ku- One may explicitly 
check that these states fulfill (II) in (15.61) . In order to fulfill (III) we find the following conditions: 



D = 2L, even 



! —X 



2L-1 

n <ik 

k=l 



X 2L + b [ M +1) ■ A 



2(L-M) 



M=l 



D = 2L + 1, odd 



(5.9) 



xX 

ru 

k=l 



\ 2L + b M L+2) ■ a2(l ~ m) 



M=l 



That is in order for (j5.9p to be satisfied, A has to be an eigenvalue, that is a root of the (real) 
characteristic polynomial 



vr 2L (A) = X 2L +J2 * } A 2(i " M) 

M=l 

vr 2L+ i(A) = X\X 2L +J2€ 0) ^ L - M) 



D = 2L even 



k = 2L + l = D + l 



(5.10) 



M=l 



D = 2L + 1 odd k = 2L + 2 = D + l 



with coefficients defined as in ( 15. 8p . Since the square- bracketed terms in the characteristic 
polynomials (15. lOf) contain only even powers of A we may replace X 2k = A k , with A = |A| 2 in order 
to arrive at a reduced purely real notation for them: 



= A 2L + C o) A 2(i " M) 



M=l 
L 



(5.11) 



M=l 



where we have pulled out a pref actor (— 1) M coming from the purely imaginary nature of the matrix 
elements (15.41) . This confirms the results obtained in [46J. 



5.4 Lower Bound on the Spectrum 

Under fairly general assumptions on the scaling of the smallest three spins with respect to the fourth 
largest spin at a vertex, an explicit lower bound for the smallest non-zero eigenvalue can be obtained 
analytically, as well as the scaling of this lower bound with respect to the maximal spin. The idea is 
to use Gershgorin's theorem [5H], which provides an upper bound for the modulus of the maximum 
eigenvalue for a general matrix in terms of the maximal row/column sums of the modulus of matrix 
elements. Now an upper bound for the eigenvalues of the inverse matrix Q^ 1 certainly provides 
a lower bound for the smallest non-zero eigenvalue of Qo- However, due to the zero contained as 
a single eigenvalue in the spectrum of Qd in the case of odd dimension D this problem, although 
presented already in [5T] . was left open there. Using certain properties of the resolvent of Qd [M], it 
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is possible to overcome this problem. For any D the spectrum of Qd can be related to the spectrum 
of the (D — l)-dimensional matrix Qd-i which is obtained by deleting the last row and column of 
the matrix in (I5.3p . As can be shown [52] the spectra of Qd and Qd-i interlace, that is between 
any two eigenvalues of Qd there is exactly one eigenvalue of Qd-i- hi case of odd dimension D, 

Qd-i is even dimensional and can be inverted without any problems, thanks to the fact that even 
dimensional matrices of the form (15.31) with non-vanishing entries on the first off-diagonal have no 
zero in their spectrum |51j . After an admittedly lengthy derivation presented in the companion paper 
[52] . a lower bound for the smallest non-zero eigenvalue ^|>4- vertex °^ ^-21) in leading order of the 
maximal spin j max is obtained as 

A^ ertex (j max ) > £W\Z\-W(l23)\- Jmax (5.12) 

Here a(123) = 0, ±2, ±4 is defined as in (13.101) . (15.21) . This may be compared to the numerical result 
Xmin = £p ■ \J\Z ■ a(123)| ■ 2 ■ a/ jmax {jmax + lj obtained in [HI] contributed by the spin configuration 

ji — 3i — \ an d j3 = 3a — jmax- Remarkably the estimate presented in [52] even explains this 
combination of spins, that is in the upper bound approximation this particular combination of spins is 
enforced in order to maximize the row/column sum of the inverse matrix and can thus be understood 
analytically. 



6 Summary and Outlook 

In this paper we present the results of an analysis of the operator corresponding to the classical 
expression for the volume of a spatial region (as constructed by Ashtekar and Lewandowski) within 
Loop Quantum Gravity. As the action of the volume operator can be decomposed into the action on 
single vertices, we can restrict our attention to a single vertex only. Consequently we have studied 
the volume operator at gauge invariant 4,5,6,7-valent vertices. The companion paper [52] contains 
further details. 

We have started from [51J, where the matrix elements of the volume operator were explicitly 
derived within a framework using recoupling theory of angular momenta. This results in a drastically 
simplified matrix expression, which is a crucial ingredient and starting point for the work presented 
here. 

For the gauge invariant 4- vertex we obtain consistency with the results obtained in [51] . Moreover 
we find an analytic derivation for a lower bound of the non-zero eigenvalues and their scaling behavior 
with respect to the maximum spin incident at the vertex. We are able to give a simple proof for 
the spectral simplicity of the single matrix expressions of the volume operator at the 4-vertex. We 
also find an expression for the eigenstates in terms of the matrix elements and the eigenvalues of the 
volume operator compatible with [IS]- In order to achieve these results we heavily take advantage 
of the band-structure of the matrix, which enables us to use certain techniques for Jacobi-matrices 
[54] . However, this implies that it will be difficult to extend these analytical investigations to the 
general matrices arising from higher valent vertices. Consequently a numerical study appears to be 
the viable way to analyze the properties of the volume operator at valences greater than 4. 

There are numerous technical difficulties to be overcome in order to do numerical computations 
for higher valent vertices. We have shown how the gauge invariant basis states can be implemented 
on a computer using the rules for the recoupling of angular momenta. A more difficult task is then 
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to understand the sign-factor combinatorics of possible embeddings of the edges of the underlying 
graph at a vertex, that is the possible relative orientations of the edge tangents at the vertex. We 
have used a Poisson sprinkling process in the actual computation, excluding vertices with coplanar 
edges. However we have also pointed towards an analytical^! solution of this problem in terms of 
coupled systems of inequalities. We could then numerically compute the matrix representations of 
the volume operator for all sets of edge spins ji < ■ ■ ■ < ]n = jmax (N = 5, 6, 7). Here we can take 
Jmax up to a certain cutoff which is determined by the limitations of the computational resources. 
The eigenvalues of these matrices are computed using singular value decomposition. The eigenvalues 
are finally collected into histograms, whose bin occupation numbers serve as a notion of spectral 
density. We finally discuss the properties of this spectral density, where we superpose the eigenvalues 
coming from all spin configurations and/or sign-factor combinations. 

Our numerical results reveal that the sign factors contained in the definition (13. 4p of the Ashtekar- 
Lewandowski volume operator open a fascinating interplay between spin quantum numbers and gauge 
invariance on one side and the graph geometry on the other side. This manifests itself in the fact 
that the spectral density of the volume operator depends on the sign configuration at a vertex, which 
is induced by the vertex geometry. This aspect has not been analyzed in detail before due to the 
already very complicated recoupling computation. 

Having obtained these results the following tasks and questions occur as important for future 
investigation. 

We have seen that the volume operator is sensitive to the vertex geometry as manifested in the 
relative orientation of edge-tangents at a vertex which result in a particular spatial diffeomorphism 
invariant sign-factor combination. Here it is important to understand what kind embeddings of N 
edges in three dimensional Riemannian space, labelled by sign factors e, can be realized at all. Here 
we also have to understand how to factor out permutations of the edges. These issues will be worked 
out in detail in a forthcoming publication [53]. It will be interesting to see how the number of realized 
embeddings of iV edges grows as N is increased, in particular if there is a kind of saturation process 
of that number as we send N to infinity, (e.g. there is the possibility that the growth in the number of 
possible embeddings rapidly decreases as the number of edges is increased.) This could also provide 
us with an answer to the question of whether it is justified to restrict oneself to graphs which have 
vertices only up to a certain maximum valence. 

In the analysis above, it appears that e.g. a 7-valent vertex can contribute relatively large volume 
eigenvalues, produced by relatively small spins, as can be seen from figure [71 Could it be a general 
feature that very high valent vertices produce very large (even macroscopic) volume eigenvalues from 
small spins, or is there some mechanism which suppresses such contributions from large valence 
vertices? 

Understanding the sign-factor combinatorics and their encoded geometrical information is im- 
portant for a number of reasons. Eventually one would like to find a combinatorial formulation of 
Loop Quantum Gravity [21]. Are the sign factors related to local geometrical quantities of cell- 
decompositions of the three dimensional spatial foliation surfaces? 

Another consideration is that it is desirable to gain more control, that is a geometric interpre- 
tation, on the quantum numbers of states in the (spatial) diffeomorphism invariant Hilbert space 
[1] of the theory. One seeks states which are, roughly speaking, constructed by an averaging of a 

31 Of course, the solutions of the inequality systems as presented in [52] use a computer algebra program. However 
the actual computations will be exact in the sense that there are no numerical errors present when working with these 
solution tables of systems of inequalities. 
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graph over all its possible (diffeomorphic) embeddings into the spatial foliation hypersurfaces. These 
diffeomorphism invariant states could be classified (in addition to the moduli parameters of [27]) by 
the volume spectrunj^l caused by the sign configurations at the vertices which are induced by the 
diffeomorphism invariant structure of the graph. 

Given the fact that we now have at hand a numerical treatment of the volume operator in full 
Loop Quantum Gravity, it will be worthwhile to look again at the suggestion pointed out in 
to use the volume operator in the construction of complexifier coherent states rather than the area 
operator. This may remove certain embedding ambiguities from which the current complexifier co- 
herent states suffer [H] , which result from the usage of the area operator as complexifier J 33 l There are 
certain conceptual issues still to be solved, however, e.g. that due to the huge kernel of the volume 
operator these coherent states are not normalizable in the construction proposed so far. Nevertheless 
this can be a promising starting point to (at least numerically) construct (spatial) diffeomorphism 
invariant coherent states for concrete computations, using the sign factors as labels to denote the 
peakedness of the thus constructed states around a given (classical) geometry. Such a notation has 
not been available up to now, which has prevented the construction of these states. 

We would like to mention here that the question of what category of diffeomorphisms and edges 
should be considered as admissible for LQG is a lively research area. In [2H] the use of the smooth 
category is discussed, while [28] looks at the possibility of using the piecewise linear category. In [30] 
diffeomorphisms which are analytic everywhere except for a finite number of points are analyzed. 
Finally, as pointed out in [31] and applied in [32J, it appears that one has to stick to the stratified 
analytic category, which roughly allows that an edge can get a finite number of kinks under such 
transformations. See also [33] on this subject. Of course, the question of what information should 
enter diffeomorphism invariant states has to be considered within this context, and it could for in- 
stance be that one has to use states which are additionally averaged over all a-configurations, and 
thus do only contain information on the linking structure of the underlying graplQ 

Another potentially important direction we would like to point at is the possible application of 
techniques from random matrix theory [55] to Loop Quantum Gravity. By inspection of the volume 
spectrum averaged over the a-configurations and spins as in figure [3J, we can see a certain similarity 
of the produced spectral densities to the densities occurring in the context of certain ensembles of 
matrices having entries which are independent random variables. The idea here is the following. The 
averaging of the volume spectrum over all spin configurations ji, ■ ■ ■ ,Jn of the iV-valent vertex v 
(up to a certain maximum spin j ma x ) can equivalently be seen as probing a set of random numbers 
jx, . . . , ]n all occurring with the same probability. The problem to be solved in this direction is to 
extend the notion of random matrices to matrices whose entries are not random variables themselves 
but functions of a set of random variables, that is the spins ji, . . . For such an analytical treat- 
ment it will be necessary to further investigate universal properties of the spectral density of the 

32 Note that the region whose volume is measured also has to be specified in a diffeomorphism invariant manner, for 
instance by a matter field. 

A surface has to be chosen which is punctured by the underlying graph, giving rise to eigenvalues of the area 
operator. This surface can be tilted with respect to the graph, such that the number of transversely intersecting edges 
changes, which makes the complexifier eigenvalues a matter of choice. 

34 From this perspective it might be interesting to look at the question of whether the operator \/(- RS ') of [33] can be 
related to the averaged V^ AL ^ analyzed in the present paper. 
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volume 



Finally we would like to point out that V is not only essential in the construction of the op- 
erator corresponding to the Hamilton constraint [7j (as shown in (12. 5ft ). but in all kinds of matter 
Hamiltonians [8] in the theory. Having achieved a computational treatment of the volume operator 
now opens the door to numerically test the results of symmetry reduced models [201 12T1 122] within 
full Loop Quantum Gravity. As the code of the present analysis will be publicly released, we hope 
that this might serve as a starting point to developing a computational toolkit for numerical Loop 
Quantum Gravity, which is freely available to the research community. 

Concluding, we would like to state that despite the conceptual problems which are inherent to 
Loop Quantum Gravity, it is important to study the physical implications of the theory in its current 
form. Numerical analysis now opens the door to new insights into the predictions of Loop Quantum 
Gravity, as it allows the study of aspects of the theory which are inaccessible using analytical tech- 
niques. The insights thus obtained provide important feedback to the development of the theory. 
To facilitate future numerical analysis, we intend to release the source code we have used in this 
investigation, as 'thorns' for the Cactus high performance computing framework. We hope this will 
a first step in the development of a community code base for a wide range of computations in Loop 
Quantum Gravity. 
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